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We consider RF-spectroscopy of ultracold Fermi gases by exact simulations of the many-body 
state and the coherent dynamics in one dimension. Deviations from the linear response sum rule 
result are found to suppress the pairing contribution to the RF line shifts. We compare the coherent 
rotation and quasiparticle descriptions of RF-spectroscopy which are analogous to NMR experiments 
in superfmid 3 He and tunneling in solids, respectively. We suggest that RF-spectroscopy in ultracold 
gases provides an interesting crossover between these descriptions that could be used for studying 
decoherence in quantum measurement, in the context of many-body quantum states. 
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INTRODUCTION 



II. COHERENT AND QUASIPARTICLE 
PICTURES 



The excitation spectrum of a quantum state can be 
measured by outcoupling particles and creating quasi- 
particles. In ultracold Fermi gases, RF-spectroscopy 
[E H, 0, 0, 0, IS 0, El] realizes such outcoupling. One 
of the two interacting (pseudo)spin components is cou- 
pled by an RF field to a third component. This can be 
viewed as a quasiparticle creation process 0, [T3, El, ED, 
0, 0, 0, EB, [13, El] analogous to tunneling. It can 
also be viewed as a coherent rotation in (pseudo)spin 
space (19L I20I l2ll. I22II . a nalogous to NMR experiments in 
superfluid 3 He [23l. 1 24 l25j |. The picture which applies 
depends on the coherence of the time evolution [22]]. We 
compare the coherent rotation and quasiparticle descrip- 
tions, and study nonlinear effects, by exact simulations 
of the dynamics in ID, using the matrix product state 
(MPS) method (H]. 



In this letter, we suggest that the amount of deco- 
herence, not the final state interactions directly, is the 
crucial issue determining which kind of theoretical de- 
scription for RF spectroscopy should be applied. We dis- 
cuss the difference between the coherent and quasiparti- 
cle pictures from this viewpoint, and make a connection 
between RF-spectroscopy and the quantum measurement 
problem [13, S|, with discussion on possible experiments. 
We give exact numerical results on how deviations from 
the linear response lead to a reduced pairing contribution 
to the RF line shifts. These predictions are directly rel- 
evant for coming ID experiments and should be a good 
estimate for the coherent 3D case as well. 



The system considered is described by the Hamiltonian 
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where U\i and Uif are the interaction strengths {U^f in- 
teraction is neglected), /ii = fj,2 = M the chemical poten- 
tials, m is the mass, rij(r) = \I/|(r)\I>j(r) and f2 the field 
(effective) Rabi frequency. For a non-interacting system 
(U12 = Uif = 0), the spectrum of transferring particles 
to the (initially empty) final state / has a peak at fre- 
quency ujfree- For an interacting system, the maximum 
transfer of atoms from state 2 to state / takes place at a 
certain RF frequency ujrf that can differ from ujf ree by 
some detuning 5 = lorf — ujf ree . 

In the quasiparticle picture, assuming BCS state with 
excitation energy E%, gap A and linear response for the 
field, the detuning is bound by the energy con servation 
relation @| 6 = E k + e k - [i = yj (e k - fi) 2 + A 2 + e k - /-*, 
which gives the threshold detuning 5th — vV + ^ — 
fx ~ A_ ~ ^A_ (minimum 5th is obtained for the non- 
interacting case with single particle energy e k — 0). In- 
cluding the Hartree interactions gives 
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where Ef is the Fermi energy and n a denotes the den- 
sities for different components. For the BCS state the 
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threshold 5th corresponds well to the maximum peak po- 
sition. Eq.(2) gives also the molecular binding energy 
5 t h ~ 2\fj,\, for n < and >> A. 

In the coherent rotation picture, within the linear re- 
sponse regime, the mean value for the frequency shift, 5, 
can be obtained from the sum rules [13, El 
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which in the BCS limit becomes |2( 
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The Hartree term (Uif — Ui2)ni is the same in Eq.(2) 
and in Eq.(4). Neglecting it one obtains for U\j = that 
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Both the quasiparticle and coherent rotation pictures 
predict that the frequency shift is proportional to A 2 , 
but for different reasons. For the coherent rotation, the 

A 2 

term is the difference between the total energy of 

the final and initial many-body states, and dividing it by 
n% gives the energy difference per particle. In the quasi- 
particle picture, the suppression of pairing by the Fermi 
energy, 5th — appears because the particles that can 
be excited with the smallest energy are at the bottom of 
the Fermi sea (ek — 0). Note that in some other contexts 
(e.g. tunneling in metals), the states at the bottom of the 
Fermi sea are Pauli blocked, whereas here the final state 
is initially empty. In fact, by having initial occupation 
of the final state, the dependence of 5th on A could be 
brought towards 5th ~ A [| H§] in the quasiparticle 
picture, whereas in the coherent rotation picture the A 2 
dependence would be unchanged. For the coherent rota- 
tion, the term (U12 — Uif) multiplying A 2 leads to the 
zero shift 5 = if the interaction strengths are the same, 
i.e. the interaction is SU(2) invariant in spin space [25| . 
In contrast, in the simplest quasiparticle picture (Eq.(2)) 
the term proportional to A 2 does not depend on the final 
state interactions. This can be valid for many-body pair- 
ing with weak interactions: a single outcoupled particle 
does not form a Cooper pair because there is no Fermi 
sea for the particles in the final state. For strong final 
state interactions and/or if bound states exists, pairing 
can take place without the Fermi sea. Such a 1-f bound 
state can be added to the quasiparticle picture to describe 
bound-bound transitions [13, LL8| . 

The shape of the spectra in the coherent rotation and 
quasiparticle pictures allows distinguishing between them 
in the experiment. In the quasiparticle picture, the spec- 
tral peak is asymmetric with a long tail. In the coherent 
rotation description, the spectra are symmetric and not 



broadened by interactions [20l . l22l |. We will now com- 
pare the coherent and quasiparticle pictures in a one- 
dimensional system, paying particular attention to the 
validity of the linear response. All the previous works in 
the literature on RF-spectroscopy have assumed linear 
response (only in [3(3] an approximative non-linear model 
was considered). We thus present the first results with 
no approximations in the field-matter interaction. 



III. METHODS AND RESULTS 

We consider the Hamiltonian (1) within the one- 
dimensional, one-band Hubbard model with nearest 
neighbour hopping and on-site interaction, assuming the 
tight binding approximation: H = —Jj2(ij)a c ia c 3<? + 

^12 E* 44^2^1 + ^1/ Ei CiAfdfCil + E i(T Vo J ■ 



Ei f (4 



Ci2 



4 f c if 



n E* (4 Ci2 + 4 c »/) i where 



c ia creates a particle a to site i. No renormalization 
is needed for the interaction strengths due to the natural 
cutoff by the first band. We set the hopping J = 1 in 
calculations, all other quantities are in units of J. The 
relation between J, U and the parameters for atoms in 
an optical lattice are as in j3l|]- The chemical potentials 
(x a are chosen so that the density becomes corresponding 
to the half filling. The final state is initially empty. 

For calculating the ground state and simulating the 
dynamics we use the matrix state product method as de- 
scribed in (26l . I32I . This formalism is especially suited for 
simulating time evolution in sufficiently regular systems. 
All interaction strengths, i.e. the BCS-BEC crossover, 
can be considered, provided the single band approxima- 
tion stays valid. The ground state (when fl — 5 = 0) is 
obtained by operating on the initial state with imaginary 
time evolution operator which we approximate by the 
Forest-Ruth formula (33|. After calculating the ground 
state, we simulate the dynamics by applying a square 
pulse. For each detuning, we calculate the expectation 
value for the number of atoms in the final state after hav- 
ing applied the pulse and thus obtain the spectrum. The 
numerically obtained RF line shift is denoted as 5 num . 

The errors due to the numerical method are controlled 
and very small: we have varied the Schmidt number in 
the MPS method from 5 to 12, and the results start to 
converge around 5 (depending on the interaction). The 
error still left in the results fits within the data point 
symbols. The simulations are double-checked using two 
computer codes (one in C++, the other in Fortran). 

Figure 1 shows the shift 5 num as a function of the inter- 
action strength, U12, when Uif = 0. In Figure 2, we vary 
Uif 7^ and keep U12 fixed. For comparison, we show 
the shift of the peak 5 as given by the sum rule within 
the linear response theory, Eq. (3), where we calculate 
(Jdr*|(r)*J(r)* 2 (r)*i(r)) numerically. Having the 
BCS result Eq. (4) in mind, we also plot the Hartree mean 
field contribution U12TI1, and extract the quantity A from 
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Figure 1: The frequency shifts S nvm for different interaction 
strengths U12 < (Uif = 0). The circles represent 1.0%, 
crosses 10% and open squares 50% particle transfer, also in 
the inset. The linear response - sum rule result 5 (solid line) 
works well for weak interactions or if the amount of trans- 
ferred particles is small. The result given by Fermi golden 
rule (filled squares) and the Hartree term Ui2Ui (dotted line) 
are also presented. The inset shows the differences between 
sum rule and numerical results 5 u num • 



Eqs. (3) and (4) via U (c^cl^cn) w Unw?, + A 2 U~ 1 . 
Using A obtained in this way we apply the standard BCS 
description to present results given by the simplest quasi- 
particle picture (Fermi golden rule). 

As Figure 1 shows, the numerical result deviates from 
the linear response sum rule result (Eq. (3)) by an in- 
creasing amount for increasing interactions, if the num- 
ber of particles transferred is e.g. 50 %, or even 10 %, 
of the particle number in state 2. Only when as little as 
~ 1 % of the atoms are transferred, the linear response 
is strictly valid. For U12 > 0, the deviation is found to 
be similiar. Figure 2 shows the effect of nonzero Uif, 
i.e. the simple amplification by the factor (Uif — U12) in 
Eq.(3). As expected for a coherent process, the spectra 
are symmetric and the width depends only on O as seen 
in Figure 3. 

We checked whether the simulated time evolution cor- 
responds to simple Rabi oscillations with some effective 
detuning caused by interactions: we fitted to the numeri- 
cal oscillations the maximum peak height s ,2+ 4 q2 and the 
period 2n(8' 2 +4f2 2 )~ 1 / 2 with 5' as the fitting parameter. 
The fitting was found to fail hand in hand with the failure 
of linear response approximation for increasing interac- 
tions. This is in accordance with the observation (see 
Figure 1) that nonlinear effects show up when the pair- 
ing contribution starts to become considerable compared 
to the Hartree term (for Hartree shifts alone, the simple 
Rabi picture should be valid). For 50 % of transferred 
particles, the suppression of the remaining pairing con- 
tribution (i.e. other than Hartree) by nonlinearity ranges 
from 30 % to 100 % for the values used in Figure 1. In 
the simplest quasiparticle picture, the nonlinear shift is 
not necessarily this dramatic since higher order processes 
might produce primarily subthreshold weight rather than 
shift of the peak. 




Figure 2: The frequency shifts S num for different final state 
interaction strengths Uif (U12 = —1.0). The circles represent 
1.0 %, crosses 10 % and open squares 50 % particle transfer, 
also in the inset. The sum rule result S is shown with the 
solid line and the Fermi golden rule result with filled squares 
(at Uif = —1 the former is and latter 0.01). The inset 
shows the differences between sum rule and numerical results 

Snum- 



The results shown here are for densities close to half 
filling. For U12 — —1, we have done simulations also for 
low densities (filling 0.25) which correspond to the con- 
tinuum ID case and found that the nonlinear deviation 
for 50% particle transfer was ~ 10% of the total shift 
(c.f. 2% for half filling and U12 — —1). We have mainly 
used 20 lattice sites but also bigger lattices give the same 
results, see Figure 3. 



IV. QUANTUM MEASUREMENT AND 
EXPERIMENTS 

It is now interesting to compare the coherent and 
quasiparticle pictures from the point of view of a quan- 
tum measurement. In the former, the system decoheres 
only after the pulse. In the latter, the measurement 
instantaneously projects the system into the total fi- 
nal state which can be e.g. one free particle in state 
/ and one quasiparticle in the superfluid. The start- 
ing point for both descriptions (within the linear re- 
sponse) is the imaginary part of the correlation function 
D f2 = -ifdr (T[v4(r,t)V/(r,t)^ (0,0)^(0,0)]) which, 
as well as the sum rules, corresponds to the coherent 
Hamiltonian evolution. On the other hand, Fermi golden 
rule corresponds to the instantaneous projection (Fermi 
golden rule result is often obtained in an experiment for a 
long, weak excitation but formally it describes an instan- 
taneous projection). Fermi golden rule follows by divid- 
ing Df2 to the final and initial state contributions. The 
measurement done is thus defined by the way the corre- 
lation function is calculated. It is possible to calculate 
Df2 by including final state interactions in a fully self- 
consistent way [2fjJ, corresponding to coherent rotation 
of all particles. In contrast, the simplest quasiparticle 
picture (Eq.(2)) is obtained by a straightforward factor- 
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Figure 3: Spectra for different values of the coupling Q, the 
interaction strength U12 and the number of lattice sites Nl. 
The shape of Fermi golden rule result is also shown (in arbi- 
trary units). The width of the spectrum depends only on SI. 
Here U\f = and the duration of the pulse is 0.8£2 _1 . 



of both quasiparticle and coherent rotation type results 
for ultracold Fermi gases, and the possibility of a tunable 
crossover between them suggest that RF spectroscopy in 
these systems may be suited for studying the long stand- 
ing problem of quantum measurement and the question 
of how decoherence sets in in the measurement process, 
in the context of many-body states. Our ID results give 
predictions for the fully coherent and the fully decoher- 
ent (quasiparticle/Fermi golden rule) values of the RF 
line shift. Another experimental indicator of the amount 
of decoherence is the asymmetry of the peak. One could 
investigate how these two observable quantities (peak po- 
sition and asymmetry) evolve - smoothly or abruptly - 
when the amount of decoherence is changed, either by 
the length of the pulse or by the strength of the final 
state interactions. One could further repeat such exper- 
iments throughout the BCS-BEC crossover to find out 
whether many-body and single particle paired states are 
vulnerable to decoherence in a different way. 



ization of £>/2- Then the spectrum obtained no longer 
satisfies the f-sum rule, i.e. it does not describe coherent 
evolution but a projection measurement to the chosen fi- 
nal state. If the final state is likely to have bound states 
or strong interactions, they should be included in the 
quasiparticle description [la. [lgl] . 

Some of the spectra obtained in the early RF experi- 
ments on pairing @, [f| 0] are symmetric although this 
could be partly due to broadening by inhomogeneous 
density in 0, @|. However, the spectral shapes depend 
on the interaction strength, density and temperature, 
and the quasiparticle picture can be applied to describe 
qualitative features such as the double peak structure 
EH, G3] ■ Exact modelling of the intermediate cases 
between the coherent and simple quasiparticle pictures - 
or realizing the experiment in such a way that one pic- 
ture is strictly valid - becomes crucial in high precision 
measurements of the pairing gap. The recently obtained 
spectra @ in the case where final state interactions have 
been reduced show clearly asymmetric shape. 

For strong initial state and weak final state interac- 
tions the eigenstate wavefunctions of the initial and fi- 
nal states are very different: the RF field then couples 
one state to a continuum of states. Decoherence is the 
strongest when a state is coupled to a large number of 
other states, which enforces a projection measurement. 
This is obviously the case for the molecule continuum 
dissociation process. If the final state resembles the ini- 
tial state, the process is more likely to stay coherent. 
This can happen e.g. for stronger final state interactions. 
The extreme case is the Hartree mean fields which do 
not affect the wave function in a homogeneous system; 
fully coherent rotation was observed in measurements of 
the Hartree shifts 0, [3| . The experimental observation 



V. CONCLUSIONS 

We have compared the exact coherent dynamics of RF 
spectroscopy of interacting Fermions in one dimension to 
the quasiparticle picture within the BCS formalism. We 
have also pointed out a potential connection with quan- 
tum measurement and decoherence. Large suppression 
of the pairing shift from the linear response prediction 
was found when the number of transferred particles was 
more than ~ 1 %, The main difference between the ID 
and 3D physics is the absence of long range order in ID. 
Since RF-spectroscopy does not measure the long range 
coherence of the many-body state, our results should give 
a good estimate of the relative effect of nonlinearity also 
in the 3D coherent case. Therefore, in 3D experiments 
where more than 1% of particles were transferred and fi- 
nal state effects were considerable (nearly coherent case) , 
it is likely that the RF shifts have been suppressed by 
amounts similar to the ones presented here. In the fu- 
ture, it would be especially interesting to consider ID 
systems where exact theory and modeling is possible, to 
understand both coherent and incoherent processes in RF 
spectroscopy. The deviations from the linear response for 
strong interactions found here also indicate that such sys- 
tems could display other interesting nonlinear phenom- 
ena in the field-matter interaction. 
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